Preprint MKPH-T- 10-02 

A dispersion relation for the pion-mass dependence of hadron properties 



Tim Ledwig, Vladimir Pascalutsa, Marc Vanderhaeghen 

Institutfur Kemphysik, Johannes Gutenberg Universitdt, Mainz D-55099, Germany 



Abstract 

We present a dispersion relation in the pion-mass squared, which static quantities (nucleon mass, magnetic moment, 
etc.) obey under the assumption of analyticity in the entire complex plane modulo a cut at negative associated 
with pion production. The relation is verified here in a number of examples of nucleon and A-isobar properties 
computed in chiral perturbation theory up to order . We outline a method to obtain relations for other mass- 
dependencies, and illustrate it on a two-loop example. 
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1. Introduction 



Present lattice QCD calculations are still limited to larger than physical values of light quark masses, iriq > iriu^d - 
5-10 MeV, but the chiral perturbation theory OifPT) [1,2] can, in many cases, be applied to bridge the gap between 
the lattice and the real world (see e.g., [3-8]). can predict at least some of the 'non-analytic' dependencies of 
static quantities (masses, magnetic moments, etc.) on pion-mass squared, or the quark mass (;«^ ~ m^). The rest of 
the contributions contain the a priory unknown low-energy constants (LECs). In this paper we examine the origins of 
non-analytic dependencies arising in ;i'PT, by considering analytic properties of the chiral expansion in the complex 
ml plane. 

The basic observation is that chiral loops exhibit a cut along the h — 

negative axis. The cut is associated with pion production which can /'' \^ 

occur without any excess of energy for < 0. Assuming analyticity in 
the rest of the m^-plane (see Fig. 1), one arrives at a dispersion relation 

of the type: i. . 

n J t -m^ + iO+ 

where / is a static quantity, 0^ is an infinitesimally small positive num- "'' 

ber In what follows, we explicitly verify this type of dispersion relation 

on a few examples of the nucleon and A(1232)-isobar properties and ^- ™t and the contour in the complex 

. , , . . t = m- plane, which go into the derivation of the 

discuss Its field of application. In particular, we consider a two-loop ex- dispersion relation in Eq (1) 
ample (a sunset graph) for which the absorptive part can relatively easy 
be extracted. We conclude by comparing this dispersion relation with a 
similar "mass-dispersion" relation long-known in the literature. 
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Figure 2: Graphs representing chiral-loop corrections to the nucleon mass. Nucleon (pion) propagators are denoted by solid (dashed) lines. 



2. Nucleon mass 

We begin right away by considering the nucleon properties as a function of f = m^. For example, the pion-mass 
dependence of the nucleon mass, computed to the nth order in the chkal expansion, can be written as: 



Ma 



= ^flffJ +^2«(f), (2) 



where a's are some linear combinations of LECs, 2^)^ (0 is the fth order nucleon self-energy given by the graphs of 
the type shown in Fig. 2. According to the power counting rules [9], a graph with L loops, A^;^ pion and A^^ nucleon 
lines, Vk vertices from the Lagrangian of order k, contributes at order p", with p being the generic Ught scale and 

„ = ^ /ty^ + 4L - 2N,: - Nn ■ (3) 

k 

The leading order pion-nucleon Lagrangian is of order k - 1, and, to the first order in the pion-field jf(x) (with index 
a - 1,2, 3), is written as [10]: 

o 

-CnN = N(x) Mn +^ {97f)T"y5]N{x) + c.t. + 0(n\ (4) 

o o ° 

where A^(jic) is the isospin-doublet nucleon field, t" are Pauli matrices, Mn, Sa, and are respectively: the nucleon 
mass, axial-coupling and pion-decay constants, in the chiral limit {m„ — » 0); "c.t." stands for counter-term contribu- 
tions, which are required for the renormalization of the nucleon mass, field, and so on. 

The self-energy receives its leading contribution at order p^, which is given by the graph Fig. 2(a) and the following 
expression: 



3^1 . r A^k k-yy5ip-y-k-y + MN)k-yy5 



'j {2nf 

4/^ 



(5a) 

p-y=Mn 



(J<?-t + iO+)[{p - kf -Ml + iO+] 

f„ (47r)2 Jq ^ ^ ' 

+ [2x2 _^2 + x)t] + ln[x2 -h (1 - x)t - /0+]) - 3Le}, (5b) 



where t = f/M^, = -1/e- 1 H-7£-ln(47rA/MAr) exhibits the ultraviolet (UV) divergence as e = (4-d)/2 — » 0, with 
d being the number of space-time dimensions, A the scale of dimensional regularization, and ^ 0.5772 the Euler's 
constant. Note that for simplicity we assume the physical values for the parameters: - 939 MeV, gA - 1.267, 
fjr - 92.4 MeV; the difference with the chiral-limit values leads to higher order eff'ects. 
After integration over the Feynman-parameter x, this result can be written as: 



with 



t^(t) = -^TTT^ir'" arccos(lV^) + Inr). (6b) 



Figure 3: Graphs representing order-p^ corrections to the electromagnetic interaction of the nucleon. 



The term in figure brackets, containing the UV-divergence, must be entirely canceled by the counter-term contribu- 
tion [11], which, to this order, is of the form: oq + 02 t, where a's contain the "bare" values of the LECs. The first term 
in brackets can be viewed as a renormalization of the nucleon mass, while the second as a renormalization of the nN 
sigma term. The remaining part, 1,^, is UV-finite and consistent with the power counting in the sense that its size is 
indeed of order p^. 

Let us now see whether this contribution obeys the dispersion relation of the type stated in Eq. (1). The imaginary 
part can be easily found from Eq. (5b) by taking into account that ln(-l + iO^) = in. 



3^: 



(47rA) 



71 

~2 



-i-t) 



3/2 



1/2 



2Mn 



(7) 



where 6 is the step function. It is quite obvious that the dispersion integral with this imaginary part diverges, which is 
consistent with the fact that the self-energy is UV divergent. From Eq. (6a) we have seen that the divergencies appear 
in the first two orders of the expansion around f = and are subsequently absorbed by the counter- terms. As the result 
one needs to make two subtractions at point f = to find 



(8) 



Now the dispersion integral, with the imaginary part given by Eq. (7), converges, and moreover, gives the result 
identical to the expression in Eq. (6b), that is, the renormalized self-energy contribution. The subtractions have played 
here the role of the counter- terms. We therefore conclude that the order-p^ self-energy correction to the nucleon mass 
obeys the suitably subtracted dispersion relation of the type of Eq. (1). We emphasize that the subtractions do not 
introduce any additional uncertainty in the result. The number of subtractions is not arbitrary but corresponds with 
the number of counter-terms available at a given order. 



3. Magnetic moment and polarizability 

We next turn to the example of chiral corrections to the nucleon's magnetic moment. For this we introduce 
the electromagnetic interaction, firstly by the minimal substitution [i.e., df^N — > df^N + 1(1 + T-ifeA^jN, d^Ti" — > 
dfji" + e"''^ eAf^n'' , with e ^ An 1131 ] in the chiral Lagrangian, and secondly by writing out the relevant non-minimal 
terms: 

= -^^^(^(1 + ^3) +K1 - ^3) k)rNF,, + c.t, (9) 

o o 

where Kp and k„ are the chiral-limit values of the proton's and neutron's anomalous magnetic moment (a.m.m.), 
respectively; furthermore y^v = tT[//Tv], F^y - d^j^Ayy According to the power-counting of Eq. (3), the leading order 
chiral correction to the electromagnetic coupling comes at order p^ and is given by the graphs in Fig. 3. These graphs 
give the following contribution to the a.m.m. of, respectively, the proton and the neutron: 



^fn)^ Jo 



(4. 
(4;rA)2 



djc- 



(1 - \x)x^ 



(\- x)(tlMl) + -iO 



Ax ■ 



x^ 



(1 -.x)(f/M2)-Hx2-/0+ 

3 



(10a) 
(10b) 



For negative t, these functions develop an imaginary part: 

lra.r-'\t) = $^2U'^r^4i'-Uir^A))ei-t), (Ua) 
lm.r''-\r) . -^.^ir.Afoi-O, (lib) 

with T = t/M^ and A = ^|t^ - r. Substituting these expression into the dispersion relation Eq. (1), and performing 
the integral, we obtain: 

II — I I / -t- T / A/ / I 

(12a) 



1 r Im<;«nO ^|A4l/, 4-llr + 3T^ ^ ^ 1 

-- d? = -< 1 Vt arccos — 6t + t(-5 + 3r)hiT> , 

-oo "y^ 4' 

1 f,^, lm^^'°°P'(0 gjMl 1 2-T ^ V? , \ 

— df = T-s2 Vt arccos rlnrV. 

nj t'-t (47rA)22\ fTTTZ 2 / 

-cx> "Y 4 



The exact same result is obtained by integrating over the Feynman-parameter in the loop expressions of Eq. (10). The 

dispersion relation proposed in Eq. (1) is thus verified in this example as well. Note that in this case we do not need a 
subtraction simply because the integral converges. However, since the complete result to this order is 

K=K+/P'^°°P^ +c.t, (13) 

the counter-term contribution, which here is just a constant involving the "bare" value of a.m.m., can be put in corre- 
spondence with one subtraction at r = 0. 

Let us remark that the same expression for the nucleon a.m.m. is obtained as well by two other dispersive methods: 
a derivative of the Gerasimov-Drell-Heam sum rule [12] and a "sideways dispersion relation" [13]. Together with the 
present result, we therefore already have three different dispersion relations, which can be applied to the a.m.m. 
calculation. One can hope that at least one of them will make the two-loop calculation of the nucleon a.m.m. more 
feasible. 

We conclude the discussion of the nucleon properties with the example of the scalar nucleon polarizabilities: 
(electric) and (magnetic). The specifics of this example is that the leading order {p^) correction comes entirely 
from chiral loops, the counter-terms are absent. In the case of magnetic polarizability of the proton, the result is given 
by [14, 15]: 

«(3).A = '^Sl r\(^ _ (l-x)(l-3x)^ + x _ xt/Ml + x^[l - (1 - x)i4 - 20x + 21x^)] | 
'^P \92n^fiMt, Jo I (1 - x){tlMl) + - iO^ [(i _ x){tlMl) + x^ - iQ^f > 

The imaginary part can be easily calculated: 

2 2 

ImySf (0 = -—fJ±—I^[2 - 72A + (418^- 246)t- (316A- 471)t2 + (54A- 212)t3 + 27r'] Oi-t), (15) 

and the dispersion relation of Eq. (1) can be shown to hold also for these expressions. The electric polarizability at 
order withstands this test too, however the expressions are more bulky and will be omitted here. 

4. A-resonance 

It is interesting to examine a case where the pion production cut extends into the physical region, as it happens in 
the case when the A(1232) is included as an exphcit degree of freedom in the chiral Lagrangian (see, e.g., [16-20]). 
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In this example the cut may extend from f - -oo up to f = A^, with A - Ma - M^, the Delta-nucleon mass difference. 
The pion-mass dispersion relation for a static quantity / becomes 



/(m^) = — df ^— — . 



(16) 



Let us demonstrate how it works on the example of a chiral correction 
to the A-isobar mass. A one-loop graph with the cut all the way up to A^ is 
shown in Fig. 4. It yields the following contribution to the self-energy [21]: 



^) f dx(xMA+MN)M\x)\L^+lnM\x)l (n) 

K I Jo L J 



y,(;rA'loop). s _ 1 / I^aMa 

where hA - 2.85 is the ttA^A axial-coupling constant, and 



M\x) 



• (1 + - t)x + r^- iO^ = (x - af 



Figure 4: A chiral loop correction to the 
A-isobar mass. Double lines denote the A 
propagation. 



(18) 



with r = Mn/Ma, t = f/M^, a - j(l + - t), and - - r^. 

The imaginary part arises again from the log, when its argument turns negative in the region of integration over 
the Feynman parameter x. More specifically. 



M\x)<0, iora-A<x<l^^^ if J <7<V . 



and as the result: 



imr 



.(;rA^ loop) 



(f) = n— 



Ma I HaMa 



\{a + r)[ - 2A' + (1 - a){T - 2A^)] + 

-^ia + r)A\ 

0, 



(19) 



t < 

< f < zf2 (20) 
t>A^. 



Note that, despite the awkward separation into regions in f, the function is continuous. For the physical value of the 
pion mass it provides the familiar expression for the A-resonance width: Fa = -2 lmZ"A - 115 MeV. 

We have checked that, analogously to the nucleon case, this chiral correction to the A mass satisfies the doubly- 
subtracted dispersion relation in the pion mass squared: 



^^(.iVloop)^ 



1 ^^^, im^r'°"'(o /»^^r 

nj t'-ml U' j ' 



(21) 



where the integration is done, as usual, in the principal-value sense. 



5. Direct calculation of the absorptive part 

Of course to find the absorptive part it should not always be necessary to go through the entire loop calculation, as 
we have done so far. In the usual dispersion relations, done in external variables such as energy, the Cutkowski rules 
offer a simple method to compute the absorptive contributions. In our case Cutkowski rules are inapplicable because 
the pion is not on it positive-energy mass shell to begin with. Nevertheless a direct computation of absorptive parts is 
possible as will demonstrated in the following three examples. 

Consider first the tadpole. Fig. 2(b), having the following generic expression 



5 



with an integer n. After the Wick rotation and going to the hyperspherical coordinates we obtain 

CO 



where J dQ4 = J^d(p j^d6 sin 6 j^dx sin^x = 2^^. The absorptive part can now be simply found as 

CO 

Im/tadCO = 7^ J dKK^^^"6iK^ + t) = T^i-t)'^" ei-t) . (24) 



(4;r)2J ' {Any 



As the second example we consider another typical integral, which appears, e.g., in the calculation of the graph 
Fig. 2(a), 

^^k 1_ 

Tf (k^-t + iO^)[(p - /t)2 - M2 + iO+] 

with = M^. Again, after the Wick rotation and the adoption of hyperspherical coordinates we obtain 



I(t,M^) = i f — ^ ^ , (25) 





and therefore the absorptive part is found as 



^WfS - ■ 



(K^ + t- jO+)(2jM cos;tr + ^ - «0+) ' 



Im/(.,M2) = -n^fdxsin'x [dK ^ s(K' + t) = L ^ ""'^ 

^' ' (2;r)4j ^ 2iMcosx + K ' {Anf ] ^ 2jM cos^ + V=i 



From these two elementary examples one can see that a direct computation of the absorptive part is in principle simpler 
than the one of the full result, since one of the integrals is always hfted by the ^-function. 

We finally come to a two-loop example, namely the pseudothreshold sunset graph with two different masses [see 
Fig. 2(c)]: 

^sunset(m^M2) =;t-'' f d'^ki f d'^ki—^ — ^ T-^ —, (28) 

j j {k]-nfi){kl-nfi)[{p-h-k2f-m 

where - M^. In this case we keep an arbitrary number of dimensions d - A - 2e, since the absorptive part has an 
ultraviolet divergence. Conveniently defining the dimensionless quantities t = in^/M^ and 

^2(2£-l) 

= HTi -Jsmsctim\ M^) , (29) 

r^(i + e) 



one can show that J satisfies a hypergeometric type of differential equation [22]: 

df2 ' L2 "'^ ' V 2' '"^'J df ' 2'^ 2e2 
As the boundary conditions one may use the easily computable massless expressions: 



- - 2. . ( - I . 4.) r ]^ 4(1 - 2e)i2 - 3e)m = ± - 2t-). (30) 



Wm^M2 = 0) = -»,2(i-2.) (31) 

e^(l - e)(l - 2e) 



J (m'-OMh - .^,2(1-2.) r(3- 4s) r(2s-i)r2(i-^)r(^) 



(32) 



hence, e.g., r^^^Jit) '^=°° -l/[s^(l - s)(l - 2e)]. 

Since for real t the equation is Hnear with real coefficients we deduce that the solution develops an imaginary part 
when the inhomogeneous term (the r.h.s.) develops an imaginary part, i.e., for t < 0. The existence of the cut for 
negative is evident. Furthermore, when seeking the solution in the e-expanded form, we again observe that solving 
for the imaginary part is simpler, because the corresponding inhomogeneous term is simpler at any given order in e, 
cf.: 

^rV-"-2) = t^-i^int+^^lnh + Oie), (33) 

ha-^r%t^-^ -2) = 6l(-0^[— -(2?- l)hi(-0 + O(e)l. (34) 
2e^ L £ J 

The solution for the imaginary part is of the form 

Im /(;) = ei-t) 7r[-^ + f(-7 + (2 + In(-O) - (1 - tf ln(l - t) + Oie)], (35) 
which agrees with the result derived from a conventional formulae [23]. 

6. Concluding remarks 

Hadron properties depend on the pion-mass squared (or, the light-quark mass, niu^d ~ ml) in an essentially non- 
analytic way. In this work we have identified the origin of this non-analyticity with a cut in the complex plane, 
which extends along the negative axis. In xP% the cut arises due to the possibihty of a "subsoft" pion production. 
Assuming analyticity in the rest of the complex plane, we are able to write a simple dispersion relation in m^, cf. 
Eq. (1). The validity of this relation has been tested here, on a number of quantities computed in;^fPT to lowest order. 
It also has been tested here on a generic two-loop example. 

There are at least two ways in which the proposed dispersion relation can be useful. First, as a consistency 
constraint of various "chiral extrapolation" formulas and methods. Second, as a computational technique, similarly to 
how the usual dispersion relations, written in terms of energy (or, in relativistic theory, the Mandelstam variables), are 
used. 

Although, the usual dispersion relation appear to be quite different from the dispersion relation in the mass, at 
a given kinematical point, the Mandelstam variables can take values given entirely in terms of mass, e.g., t = m^. 
Then, the dispersion relation in that variable can be used to connect to another kinematical point, e.g., f = 0. This is 
precisely the strategy that had long ago been proposed to relate the scattering amplitudes at the physical pion mass 
with their chiral-hmit values [24, 25]. It led to the so-called "mass-dispersion relations" [26-28], which appear to be 
similar to the relation put forward in this work.' Whether the similarities extend beyond the general form is not easy 
to tell since the analytic properties of the pion-mass dependence appear to be much more involved in the case of the 
old mass-dispersion relations. Independently of whether there is a connection between the old and new pion-mass 
dispersion relations, their test and validation in;^'PT has been made only now. 

Acknowledgments 

We thank Akaki Rusetsky for valuable remarks on the manuscript. This work is partially supported by Deutsche 
Forschungsgemeinschaft (DFG). The work of T. L. is partially supported by the Research Centre "Elementarkrafte 
und Mathematische Grundlagen" at the Johannes Gutenberg University Mainz. 



'Comparing Eq. (2.1) in Ref. [24] with Eq. (1) here, one can see that the generic form is the same. 
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